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Abstract 

This work is concerned with the formulation of the boundary quantum inverse scattering method 
for the xxz Gaudin magnet coupled to boundary impurities with arbitrary exchange constants. The 
Gaudin magnet is diagonalized by taking a quasi-classical limit of the inhomogeneous lattice. Using 
the method proposed by Babujian, the integral representation for the solution of the Knizhnik- 
Zamolodchikov equation is explictly constructed and its rational limit discussed. 
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1 Introduction 



Intcgrable quantum field theories with boundaries have been subject of intense study during the past 
decades. The great interest in such theories stems from the large number of potential applications in 
different areas in physics, including open strings, boundary conformal field theory, dissipative quantum 
phenomena and impurity problems. 

The Gaudin magnet has its origins in [1] as a quantum intcgrable model describing N spin-i particles 
with long-range interactions. The Gaudin type models have direct applications in condensed matter 
physics. They also have been used as a testing ground for ideas such as the functional Bethe ansatz (BA) 
and general procedure of separation of variables [2, 3, 4]. 

The model proposed by Gaudin was later generalized by several authors [5, 6, 7]. The spin-s XXY 
Gaudin model was solved in [8] by means of the off-shell algebraic BA. 

The XYZ Gaudin model was constructed and solved in [3] and [9] by means of the algebraic BA 
method. The boundary XXY spin-i Gaudin magnet was investigated by Hikami [10] and the Gaudin 
models based on the face-type elliptic quantum groups and boundary elliptic quantum group, as well as, 
the boundary XYZ Gaudin models were studied in [11] by means of the boundary algebraic BA method. 
In [12] the XXZ Gaudin model was solved with generic integrable boundaries specified by generic non- 
diagonal i-T-matrices. 

The Knizhnik-Zamolodchikov (KZ) equations were first proposed as a set of differential equations 
satisfied by correlation functions of the Wess-Zumino-Witten models [13]. The relations between the 
Gaudin magnets and the KZ equations has been studied in many papers [8, 14, 15, 16, 17, 18]. In [10], 
Hikami gave an integral representation for the solutions of the KZ equations by using the results of the 
boundary Gaudin model. 

In addition, the quantum impurity problem, which has been extensively investigated with renormal- 
izing group techniques [19] and conformal field theory [20], is also very interesting in itself. Andrei and 
Johannesson [21] first considered an impurity spin-s embedded in an integrable spin-i XXX chain with 
periodic boundary conditions. Subsequently, Schlottmann et al [22] generalized it to the arbitrary spin 
chain. The standard approach to dealing with the impurity integrable problem is also the algebraic BA 
method. The Hamiltonian of the impurity integrable spin chain can be constructed from the inhomoge- 
neous transfer matrix. The key point is to find some inhomogeneous vertex matrices, which satisfy the 
same Yang-Baxter relation of the homogeneous matrices, corresponding to impurity spins. 

We note that Sklyanin [7], Mezincescu and Nepomechie [24] have used a constant number K matrix 
to construct their model, where the K matrix induces the boundary fields and boundary bound states 
[23, 25]. In [26, 27], Wang and coworkers first introduced the operator K matrix to study the Kondo 
problem in one-dimensional strongly correlated electron systems. In a previous paper [28], the problem 
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of an open spin-i Heisenberg chain coupled to two spin-s impurities sited at the ends has been studied. 
Following this idea, Shu Chen at al [29] have considered de XXZ chain coupled to impurity spins with 
different coupling constants on the boundary. 

In this paper, we continue to study the XXZ chain in order to establish a link between the Gaudin 
model and the impurity problem. In the first part we construct the eigenstate of the Gaudin magnet with 
impurity by taking a quasi-classical limit of the transfer matrix for the inhomogeneous open spin chain. 
The Hamiltonian is given as a solution of the classical Yang-Baxter equation. In the second part, the KZ 
equation is studied in the context of the impurity problem; the integral solution of the KZ equation is 
obtained in terms of the Bethe eigenstate of the Gaudin magnet with impurity. 

This paper is organized as follows. In section 2 the boundary algebraic Bethe ansatz is reviewed. The 
transfer matrix for su(2) spin-i XXZ chain is construct in terms of the i?-matrix and the operator K 
matrix. In section 3 we find the Hamiltonian of the Gaudin magnet with impurity by taking a quasi- 
classical limit of the double row transfer matrix with operator K matrices and subsequently its spectrum 
is obtained. In section 4 the off-shell Bethe ansatz of Babujian is used to find the explicit integral solution 
of the KZ equation. The section 5 is reserved to a summary and discussion. 

2 The quantum inverse scattering method 

It is well know that in an integrable problem the quantum i?-matrix satisfies the Yang-Baxter equation 
(YBE): 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u - v) (2.1) 

As usual, Rij means the matrix on ® ® V^acting on the ith and jth spaces and as an identity 
on the other spaces. The variables u and v are called the spectral parameters. As a solution of (2.1), we 
use the i?-matrix for the six- vertex model defined as 

/ 1 \ 

, . _ b{u,r]) c{u,r]) 

n ^ r '>- C(U,7]) b( U ,7]) 

\ I ) 

where 



(2.2) 



sinh it sinh?7 

b ( u ^) = ^-o — ; — \ ' c ( u . f ?) = ^-r? — : — v ( 2 - 3 ) 

smn(u + T)) smn(u + rj) 

Note that in addition to the spectral parameter u, we have a deformation parameter rj which parametrizes 
the anisotropy. Moreover, this solution has the following properties: 

regularity : R(u = 0, rj) = P 
quasi — classical condition : R(u, 77 = 0) = 1 (2.4) 
unitarity : R(u, r])R(— u, rf) = 1 
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where P is the permutation operator: P \a) <g> \/3) = \(3) ® \a) and U denotes the transposition of the ith 
space. 

Let us define the monodromy matrix T(u) for the inhomogeneous TV-sites spin chain, introducing the 
inhomogeneity in the lattice through the parameter zn € C, by 

T (u\z) = R 0N (u -z N )- R 02 {u z 2 )R 01 (u - Zl ) = ( J^g ) (2.5) 

Here the operator matrix elements act on the full Hilbert space V® N . Due to the additive property of 
the spectral parameter, the YBE also holds for the inhomogeneous lattice and we have 

Riz(u - v) [T{u\z) ® T(v\z)} = [T(v\z) ® T(u\z)\ R 12 (u - v) (2.6) 

Equation (2.6) gives the fundamental algebraic structure for the QISM and gives us the commutation 
relations between the operators A(u\z),B(u\z),C(u\z) and D(u\z) 

[B(u\z),B(v\z) 

A(u\z)B{u\z) = - 1 , B{v\z)A(u\z)- C ^~ U \ b{u\z)A{v\z) 
' ' ■ b(v — u) 
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b(v- 


u) 


1 




b(u- 


v) 


c{u — 


v) 



D{,i\i)B(r\i) - 77 — —B(v\z)D(u\z) - C ^_ V J b(u\z)D{v\z) 
[B(u\z),C(v\z)} = ^—^.(D{v\z)A{u\z)-D{u\z)A{v\z)) (2.7) 
To construct the eigenstatc of our system, we can use the reference state |0) 

It is easy to check that this state is the eigenstate of the operators A(u\z) and D(u\z), and annihilated 
by C{u\z): 

A(u\z)\Q) = \Q), D{u\z) |0) = n s Jh(H-^+ ) r ? ) |0) ' C ^ z) |0) = ° (2 ' 9) 

We note that the Bethe state, a sum of spin waves [30], can be generated by operators B(u\z) acting on 
the reference state |0) . 

In order to construct an integrable open chain with boundary impurities, it is necessary to introduce 
the reflection matrices K~{u) and K + (u) which satisfy the following reflecting equations [7]: 

R 12 (u - v)K{(u)R 12 (u + v)K 2 {v) = K^{v)R 12 {u + v)K^(u)R 12 (u - v) (2.10) 

Ri2(~u + v)K+ l1 (u)R 12 (-u -v- 2rj)K+ t2 (v) = K+ t2 (v)R 12 (-u-v- 2r 1 )K+ tl (u)R 12 (-u + v) 

(2.11) 
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The solution Kf = K ± ® I and Kf = I ® K ± are the simplest reflection matrices which satisfies the 
reflecting equations. The inhomogencous transfer matrix t(u) is defined as 

t{u\z) = Tr {K+{u)T {u\z)K^{u)T^\-u\z)) (2.12) 

and forms a one-parameter commutative family 

[t(u\z),t(v\z)] = (2.13) 

where the monodromy matrix T{u\z) is given by (2.5) and by virtue of the unitarity property of our 
R- matrix (2.4), it follows the expression for the reflected monodromy matrix 

T ~ 1 (-u\z) = R i(u + z 1 )R 02 (u + z 2 ) ■ ■ ■ Ron(u + z N ) = 

It was proved in [29] that if t obeys the Yang relation 

Ri 2 (u - v)ti(u)t 2 (v) = t 2 {u)ti(v)Ri 2 (u - v) (2-15) 

then 

K~(u) =t{u + c)t- 1 (-u + c) (2.16) 

also obeys the reflecting equation (2.10) and c is a constant. 
Therefore we can construct our reflection matrix as 

Kq(u) = R 0L (u + c L )R-l(-u + c L ) = R 0L (u + c L )R 0L (u - c L ) (2.17) 

where cl is a constant decided by the left boundary. This construction give us an operator if -matrix 
instead of a constant numerical matrix, where it is identified as impurity and is not a pure reflection. 
In order to obtain the integrable Hamiltonian in the open chain we define 

U a (u\z) = K+(u)M {u\z)Kq{u)Mq\-u\z) (2.18) 

where Kq(u) = I and Kq(u) is given by (2.17), and M (u\z), Mq 1 (-u\z) are defined as 

M (u\z) =R 0R (u + c R )T (u\z), M l (u\z) = T \-u\z)R QR {u - c R ) (2.19) 

Here c R is a constant decided by the right boundary. Now one can prove that U (u) satisfies the reflecting 
equation 

R 12 (u - v)Ui(u)R 12 (u + v)U 2 {v) = U 2 (v)R 12 (u + v)Ui(u)R 12 (u - v). (2.20) 



A(u\-z) B(u\-z) 
C(u\-z) D(u\-z) 



iO 1 A\ 



4 



In order to derive the corresponding Hamiltonian we first recover the homogeneous case taking z% = 
0, i = 1, ..N and the new transfer matrix 

X(u) = Tr(K+(u)Mo(u)K (u)M^(-u)) (2.21) 
The Hamiltonian can be obtained by 

d N_1 

3 = 1 

+ coshcL + <j i <j l) + cosh 77 a\a z L 

+ coshcR {o N o\ + cr^cr^) + cosh?7 a z N a z R (2.22) 

This Hamiltonian describes, beside the bulk, the interaction of the particles with the impurities in 
the magnetic system. The contribution of the right and left impurities is explicitly in the form of the 
Hamiltonian. 

In the following we will use the algebraic Bethe ansatz [7] to solve the spectrum of X(u) for the 
inhomogeneous TV-site spin chain. 

The double row monodromy matrix can be written as 

Mu\z) = ( ) (2.23) 



C(u\z) V{u\z) 

Note that the operators A,B,C and V act on the Hilbert space V® N . The boundary Yang-Baxter relation 
(2.10) can be rewritten in terms of the operator matrix elements of U a (u\z).For convenience we define 
the operator 

V{u\z) = sinh(2w, + rf)V{u\z) - sinhr/ A(u\z). (2.24) 

then we will have the following commutation relations 

[B(u\z),B(v\z)} = (2.25) 

i nv,/ i n sinh(w + v) sinh(u — v — rf) „, 

A(u\z)B(v\z) = —4 '- -Vt7 4-B(v\z)A(u\z) 

v 1 ; v 1 ' sinh(w + v + rf) sinh(w - v) v 1 ; v 1 ; 

sinh(2w)sinh?7 

B(u|,z).4(u|z) 



sinh(2w + 77) sinh(w — w) 
sinhr/ 

sinh(u + v + 77) sinh(2i; + 77) 



£(u|z)£>(u|z) (2.26) 



1 . „ , . . sinh(u — 7; — 77) sinh(u + v + 2rf) „ . . a . , 

D Ti z B tj z = A— - - v -+-B(v\z)V(u\z) 

K 1 ; v 1 ' sinh(?j - tj) sinh(w + v + 77) v 1 ; v 1 y 

sinh(2?j + 27?) sinh?7 „. . a . , 
' -B(u\z)V(v\z) 



sinh(27; + 77) sinh(ti — v) 



S inh77sinh(2,)sinh(2, + 277) B 

sinh(u + u + 77) sinh(2i; + 77) v 1 ' y 1 ; v ; 
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The transfer matrix X(u) can be expressed as 
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and 

M 

* a (M) =B(u\z)l[B(u b \z)\0) (2.36) 

The relation (2.33) shows that the Bethe state ^({w}) is an eigenstate of the transfer matrix X(u) 
under the condition T a = 0, a = 1, M, i.e. 

-j^r sinh(u a — Zb + rj) sinh(w a + z h + rj) _ cosh 2 (u a + rj) 
1J- sinh(u a - Zb) sinh(u a + Zb) cosh 2 u a 

x rr sinh(u a + Cj) sinh(w a - Cj) yr sinh(u a - u b + rj) sinh(tt a + u b + 2r/) ^ 
j=£- R sinh(u a + Cj + rj) sinh(u a - Cj + rj) A A s inh(u a - Ub - rj) sinh(u Q + u b ) 

This equation with Zk = corresponds to the Bethe ansatz equation for XXZ spin chain with impurities 
derived in [29]. 

3 The Gaudin magnet 

We will show how the Gaudin magnet with impurity can be derived from the identity (2.33). The 
Gaudin magnet can be obtained by taking the quasi-classical limit rj — > of the transfer matrix for the 
inhomogeneous spin chain [1]. This fact indicates that the Hamiltonian is written in terms of the solution 
of the classical YBE. 

Due to the quasi-classical condition, we have the power series expansion around the point 77 = for 
each term in (2.33,2.34). 



X(u = z 3 ) = -ri + ^Hj +o(ri 3 ), A(u = Zj) = -rj + ifE, + o{r, z ), (3.1) 
and from (2.35) we have 

2 2sinh2jcoshw 3 

Fa = -V ^-T7 r^TT 1 \f a + °^ ) ( 3 ' 2 ) 

smh(zj — u ) smh(zj + u a ) 

where we have the Hamiltonian Hy. 

H > 1 = g sinh(z - + z k ) {^t + °t a > + 5 C ° Sh( ^ + Zk) + } 

W 1 f 1 

+ E sinh(z 3 - z fc ) + a t°k + 2 C ° Sh( ^ ~ Zk) + ^ 

+ ^T7 2 &1 " h Z f 7 r (coshc L (aj at + c^+cr") + ^cosh(^) (cr^crf, + 1)1 

smh(zj ■- c L )smh(x 3 + c L ) 1 ^ 1^) 2 w/ v j ^ ; J 

1 c ~T^~b i \ jcoshc^ ((7^(7+ + cr^C77) + 1 cosh(^) (cr^cr- + l) 1 , 

sinh(z 3 - cr) sinh(z 3 + cr) [ 2 J v J 7 J 

(3.3) 



+- 
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the energy Ey. 



E i = ^UO, \ + Yl l coth (zj - Cfc) + coth( Zj + c fc ] 



M 



- J2 [coth( Zj - v a ) + coth{ Zj + «„)] (3.4) 



0=1 

and the unwanted factor / a : 

f a = 2 tanh v a - ^ [coth(u a - c fc ) + coth(u a + c k )] 

k=L,R 

M N 

+2 ^2 [coth(w a - Vb) + coth(v a + v b )] - ^ [coth(w Q - z k ) + coth(w Q + z k )] (3.5) 
6=1 fe=i 

Ma 

Note that we have used the notation (u = Zj) to mean the residue at u = Zj. The integrability of these 
Hamiltonians follow from their commutativity 

[H,,H k } = 0, j = l....,N (3.6) 

which is obtained from the commutativity of the transfer matrix. 

The Bethe states "3/ (2.32) and ^ a (2.36) have the following expansions 

*(«) - V M 4> + o(r 1 M+1 ), ^ a {z 3 ) = ^-VTf + o^) (3.7) 



with 



M ( N 



and 



b^a 

where Z a are the impurity contributions 

Z a = I -r-r-r- r + ■ , , \ r r + 



sinh(w a - cl) sinh(w a + c L ) ) L Vsinh(w a - c R ) sinh(w a + c R ) J R ,g 1q n 

When we combine the terms proportional to rj M+1 in (2.33), we obtain the so-called off-shell Bethe 
ansatz equation 

1 „ , 2 sinh z,- cosh t> a ,„ . ,„,,n 

^ sinh(zj - u ) sinh(z 3 + v a ) ■? 
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This equation suggests the Bethe state </> (3.8) is an eigenstate of the Gaudin's Hamiltonian Hj, if the 
set of rapidities {v a } satisfy f a = (a = 1, M), i.e. the Bethe equations: 



N 



^ [coth(w Q - z k ) + coth(u a + z k )} 
k=i 

M 

= 2tanhu a - [coth(v a - c k ) + coth(u a + c k )] + 2 [coth(v a - v b ) + coth(u a + v b )} 

(3.12) 



k=L,R b=l 



We have derived the eigenstate and the energy of the XXZ-type Gaudin magnet coupled with impurity 
spin with different coupling constants cl and cr on the boundary. 

4 The Knizhnik-Zamolodchikov equation 

We consider the KZ equation 

V^ = 0, j = l,...,N (4.1) 
where the differential operator Vj is defined by use of Gaudin's Hamiltonians Hj (3.3): 

V, = n± - Hj (4.2) 
and k is an arbitrary parameter. The integrable condition for a set of the KZ-type differential operators 

[Vj,V k ]=Q for j,k = l,..N (4.3) 

is satisfied by the condition 

dH 1 = dJh (4 4) 

dz k dzj 

Following the idea of [8], we define the hypergeometric function A"(w|2;) by a set of differential equations 

dX(v\z) 



dzj 



= EjX(v\z), j = l,...,N 



; '' Viiy ' ] = f a X(v\z), a=l,...,M (4.5) 



dv a 

The integrability of these differential equations follows from the conditions 

dEj _ dE k dEj _ d.f a df a _ df b 



dz k dzj ' dv a dzj ' dv b dv a ^ ^ 
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In fact, it is straightforward to solve the differential equations (4.5); its solution is the function 

M 

X(v\z) = F[(coshi; a ) 2/K [sinh(w a - cl) sinh(v a + cl) sinh(w a - cr) sinh(u a + cr)]~ 1/k 

a 

iV 

x ^Q(tanhz J ) 1/ ' 2K [sinh(z J - c L ) sinh^j + c L ) sinh(zj - c R ) sinh(^j + cr)] 1/k 

3 

N M M 

x \\ - v a ) sinh(zj + v a )]- 1/K JJ[sinh(w a - v b ) sinh(w a + v b )] 2/K 

j a a<b 

(4.7) 

One can introduce the wavefunction &(z) in a integral form, which has a hypergeometric kernel [31], 

as 

- M 

*(*) = f Y[dv a X(v\z)<f>(v\z) (4.8) 

The integration path C is taken over a closed contour in the Riemann surface such that the integrand 
resumes its initial value after v a has described it. The integral function ^(z) is in fact a solution of the 
KZ equation (4.1). 

To prove (4.1) we use the fact that the Bethe state <f> (3.8) satisfies 

®?L-X^( cosh (^ Z costing + z j ) \ 

d Zj ~ ^ Uinh 2 K - Zj ) sinh 2 K + Zj )J 9a [ ™> 

where <j) a is defined in (3.9). One sees that the function 4> a does not depend on v a . Then equality (4.1) 
can be verified as 

M . , r, : \ M 



dzj 



I {€/ + «* J£) n *■ - £ + n *. 

, , / / 2 sinh Zj cosh w a _ ^,d<t> \ tt , 

= - ^ ^ sinh( , . _ Va) sinh( , . + ^ « Wa °7** - ) n *>• 

. « /" <9 / , , 2 sinh z,- cosh w a \ , , -^4 , 

= -gi ^ - J,,,,,;, + J n *» 

= (4-10) 

5 Discussion 

We have constructed and solved a Gaudin magnet with impurity. The integral representation of the 
solution for the corrcspondig KZ equation was obtained and its rational limit is here given by 
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M N 



m = /^n[^- c !)(^-4r 1/K n^) 1/2K [(^ 2 -c!)(^-4)]^ 



c 

N M M 

x nn^-^r i/K n^-^ 2/K 

j a a<b 



X 



M 

where dv — Y[ dv a . 

a 

Now, if we use our Bethe reference state (2.29) as a particular chiral primary field of the defect 
conformal field theory [32] based in the su(2) WZW conformal field theory we conjecture that the integral 
representation (5.1) should be the candidate for the corresponding M-point correlation function since it 
is a solution of the KZ equation 4.1. Here we note that the corresponding result without impurity was 
already obtained by Hikami [?]. 

In this paper we only considered the su(2) XXZ spin chain case. The generalization to the su(n) 
should be done from the view point of the Gaudin magnet with impurity [33]. Another also interesting 
case to be considered is the elliptic XYZ Gaudin magnet with impurity. 
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